Onset of Glassy Dynamics in a Two-Dimensional Electron System in Silicon 
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The time-dependent fluctuations of conductivity a fiave been studied in a two-dimensional electron 
system in low- mobility, small-size Si inversion layers. The noise power spectrum is ~ 1//°^ with 
a exhibiting a sharp jump at a certain electron density ris = rig. An enormous increase in the 
relative variance of a is observed as ris is reduced below n^, reflecting a dramatic slowing down of 
the electron dynamics. This is attributed to the freezing of the electron glass. The data strongly 
suggest that glassy dynamics persists in the metallic phase. 
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The possibility of a metal-insulator transition (MIT) 
in two dimensions (2D) has been a subject of intensive 
research in recent years [0,^ but the physics behind this 
phenomenon is still not understood. It is well established 
that the MIT occurs in the regime where both Coulomb 
(electron-electron) interactions and disorder are strong. 
Theoretically, it is well known that, in the strongly lo- 
calized limit, the competition between electron-electron 
interactions and disorder leads to glassy dynamics (elec- 
tron or Coulomb glass). Some glassy properties, such as 
slow relaxation phenomena, have been indeed observed in 
various insulating thin films |^-||]. Furthermore, recent 
work has suggested that the critical behavior near 
the 2D MIT may be dominated by the physics of the in- 
sulator, leading to the proposals that the 2D MIT can 
be described alternatively as the melting of the Wigner 
glass or the melting of the electron glass |^]. It is 
clear that understanding the nature of the insulator rep- 
resents a major open issue in this field. Here we report 
the first detailed study of glassy behavior in a 2D system 
in semiconductor heterostructures. The glass transition 
is manifested by a very abrupt onset of a specific type 
of random-looking slow dynamics, together with other 
signs of cooperativity. Our results strongly suggest that 
the glass transition occurs in the metallic phase as a pre- 
cursor to the MIT, in agreement with recent theory [ p^of . 

While glassy systems exhibit a variety of phenom- 
ena [|ll|, studies of metallic spin glasses have demon- 
strated that mesoscopic, i. e. transport noise mea- 
surements are required in order to provide definitive in- 
formation on the details of glassy ordering and dynamics. 
Fluctuations of conductivity a as a function of chemical 
potential (or gate voltage V^, which controls the carrier 
density ris) have been investigated extensively in the in- 
sulating regime [fL3[ and near the MIT [H in a 2D elec- 
tron system in mesoscopic Si metal-oxide-semiconductor 
field-effect transistors (MOSFETs). The latter study in- 
dicated that Coulomb interactions dominate the physics 
near the MIT. In order to get reasonably reproducible 
fluctuations as a function of Vg , it was necessary to make 
very slow sweeps of many hours over a very narrow range 
of Vg. Thus, all measurements represented a time aver- 



age. As a matter of fact, it had been already known jl^] 
that, at fixed Vg (or n^), a fluctuates as a function of 
time. Both high- and low-frequency fluctuations were 
evident. It was speculated that the time dependence of 
(7 was due to correlated transitions of electrons between 
different configurations jl^] or, in other words, between 
different metastable states in the glassy phase, but there 
has been no detailed study of these effects up to now. 
Here we present the first systematic study of transport 
and noise in a strongly disordered, mesoscopic 2D system 
over a wide range of and T. 

Most of the measurements were carried out on a 1 jiia 
long, 90 /im wide rectangular n-channel Si MOSFET 
with the peak mobility of only 0.06 m^/Vs at 4.2 K 
(with the applied back-gate bias of —2 V). The sam- 
ples were fabricated using standard 0.25 /im Si technol- 
ogy [|l^ with poly-Si gates, self-aligned ion-implanted 
contacts, substrate doping Na ~ 2 x lO^^cm"^, ox- 
ide charge Nqx = 1.5 x lO^^cm"^, and oxide thickness 
dox = 50 nm. The fluctuations of current I {i. e. a) were 
measured as a function of time in a two-probe configura- 
tion using an ITHACO 1211 current preamplifier and a 
PARI 24 A lock-in amplifier at ~ 13 Hz. The excitation 
voltage Vexc was kept constant and low enough (typically, 
a few /iV) to ensure that the conduction was Ohmic. A 
precision DC voltage standard (EDC MV116J) was used 
to apply Vg. The current fluctuations as low as 10~^^ A 
were measured at 0.13 < T < 0.80 K in a dilution re- 
frigerator with heavily filtered wiring. Relatively small 
fluctuations of T, V^, and Vexc were ruled out as possi- 
ble sources of the measured noise, since no correlation 
was found between them and the current fluctuations. 
In addition, a Hall bar sample from the same wafer was 
measured at T = 0.25 K in both two- and four-probe 
configurations, and it was determined that the contact 
resistances and the contact noise were negligible. 

The relative fluctuations {a — (a)) /{a) (averaging over 
time intervals of several hours) are shown in Fig. |l| for 
a few selected ris at T = 0.13 K. It is quite striking 
that, for the lowest n^, the fluctuation amplitude is of 
the order of 100 %. In addition to rapid, high-frequency 
fluctuations, both abrupt jumps and slow changes over 
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FIG. 1. Relative fluctuations of a vs. time for different Us 
at T — 0.13 K. Different traces have been shifted for clarity, 
starting with the lowest ris at bottom and the highest at top. 
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FIG. 2. (cr) vs. T for different Us. The data for many other 
Us have been omitted for clarity. The error bars show the size 
of the fluctuations, n*, n^, and Uc are marked by arrows. 
They were determined as explained in the main text. 

periods of several hours are also evident. The amplitude 
of the fluctuations decreases with increasing either or 
T, as discussed in more detail below. 

Figure || shows the time-averaged conductivity {cr) 
as a function of T for different Ug. The behavior of 
{a{ns,T)) in our samples is found to be somewhat simi- 
lar to that of high-mobility Si MOSFETs. At the high- 
est ns, for example, our devices exhibit a metalliclike 
behavior with d{a)/dT < 0. The change of {a) in a 
given T range, however, is small (only 6% for the highest 
Us = 20.2 X lO^^cm-2) as observed in other Si MOS- 
FETs with a large amount of disorder 00- d{a)/dT 
changes sign when (<j(n*)) = 0.5 e'^/h similar to other 
2D samples [|]. Even though the corresponding density 
n* = 12.9 X lO^^cm"^ is much higher due to a large 
amount of disorder in our devices, the effective Coulomb 
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FIG. 3. a) (cr) VS. T"^ for several n^lO^^cm"^) in the 
insulating regime. The error bars show the size of the fluc- 
tuations, and the lines are fits to (a) oc exp(— To/T). Inset: 
To vs. Us with a linear fit, and an arrow showing ric. b) 
(cr) vs. T^'^ for a few ris (lO^^cm"^) near Uc. The solid lines 
are fits; the dashed line is a guide to the eye, clearly showing 
insulating behavior ((cr(T ^ 0)) = 0). 



interaction is still comparable to that in other 2D systems 
{ts ~ 4, r^-ratio of Coulomb energy to Fermi energy). 

The density n*, where d{a)/dT = 0, has been usu- 
ally 0] identified with the critical density for the MIT. In 
high-mobility Si MOSFETs, the critical density has been 
also determined [|2l[ as the density ric where activation 
energy associated with the insulating exponential behav- 
ior of (cr(T)) vanishes. It was established that n* ~ nc, 
although a small but systematic difference of a few per- 
cent has been reported |ll7| , pQ[ such that n* > ric. For 
the lowest Us in our experiment, the data are best de- 
scribed by the simply activated form (a) oc exp(— Tq/T) 
[Fig. ^(a)], consistent with other studies close enough to 
the MIT [Q. The data could not be fitted satisfactorily 
to any variable-range hopping (VRH) law regardless of 
the prefactor |2^]. To decreases linearly with increasing 
Us (Fig. ||(a) inset), and vanishes at Uc ~ 5.2 x lO^^cm"^. 
Close to nc, the data are best described by the metallic 
power-law behavior {a{ns^T)) = a{ns) + h{ns)T^ with 
X ~ 1.5 [Fig. ^(b)]. The fitting parameter a{ns) is rel- 
atively small and, in fact, vanishes for ns(10^^cm~^) = 
4.72 and 4.92. Such a simple power-law T-dependence of 
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FIG. 4. Sa/{a) (main) and ^cr (inset) vs. ris at different T 
(•: 0.130 K, x: 0.196 K, □: 0.455 K, triangle: 0.805 K). Main: 
n*, rip, and Uc are marked by arrows. The vertical dashed line 
shows the region of densities ris < rig where various glassy 
properties have been observed. Inset: The vertical dashed 
line shows the location of the critical density Uc. 



cr, given by {(j{nc^ T)) oc T^, is consistent with the one ex- 
pected in the quantum critical region (QCR) of the MIT 
based on general arguments |2^], and with the behavior 
observed in 3D systems [|| and other Si MOSFETs |] 
within the QCR. Therefore, based on the analysis of 
(cr(ns,T)) in both insulating regime and QCR, we con- 
clude that the critical density ric = (5.0 ±0.3) x 10^-^cm~^ 
{ts ~ 7), which is more than a factor of two smaller 
than n*. Such a large difference between ric and n* is 
attributed to a much higher amount of disorder in our 
samples than in high-mobility Si MOSFETs 0,|o|,|l|. 

The fluctuations of a with time have been studied first 
by analyzing Sa = {{cr — (a))^)-^/^. Fig. || inset shows 
that, while da does not seem to depend on T, it increases 
with Us by three orders of magnitude. The most striking 
feature of the data, however, is the sudden and dramatic 
change in the rate of its Us dependence (see "kink" in 
Fig. Q inset), which occurs near Uq. Although we are not 
aware of any theoretical work relevant to this problem, we 
note that the observed 6a {ris) is plausible: once electrons 
enter the localized {i. e. insulating) phase by reducing 
rig below nc, their ability to change configurations will be 
severely impaired, resulting in a much more rapid drop 
of the fluctuation amplitude da with decreasing n^. 

The main part of Fig. ^ shows that, for high n^, 
the relative amplitude of fluctuations Sa/{a) is indepen- 
dent of ris and T. However, below a certain density 
rig = (7.5 ±0.3) X 10^-^cm~^, which does not seem to de- 
pend on T, an enormous increase of Sa/{a) is observed 
with decreasing either ris or T. It is interesting that 
Sa/{a) does not exhibit any special features near ric or 



n*. The onset of strong noise at rig > ric is, in fact, con- 
sistent with the observation [|2j] that, in some materials, 
a considerable increase in noise occurs in the metallic 
phase as a precursor to the MIT. We show below that 
here rig represents the density below which the 2D elec- 
tron system freezes into an electron glass. 

The normalized power spectra Si{f) = S{I^f)/l'^ (/- 
frequency) of {a — {a)) /{a) were also studied for all 
ris and T. Most of the spectra were obtained in the 
/ = (10~^ — 10~^) Hz bandwidth, where they follow the 
well-known empirical law Si = [ p5| , p6[ {p is in- 

versely proportional to the number of fluctuators and, 
usually, a ^ 1). The background noise was measured by 
setting / = for all ris and T. It was always white and 
usually several orders of magnitude smaller than the sam- 
ple noise. Nevertheless, a subtraction of the background 
spectra was always performed, and the power spectra of 
the device noise were averaged over frequency bands (^ 
an octave). Some of the resulting Sj are presented in 
Fig. ^(a). At the highest ris (not shown), Si{f) does not 
depend on ris- However, it is obvious that, by reducing 
ris below n^, Sj increases enormously, by up to six or- 
ders of magnitude at low /. This striking increase of the 
slow dynamic contribution to the conductivity is consis- 
tent with the behavior of Sa/{a) (Fig. ||). In fact, since 
relative variance {da)'^/{a)'^ = J Si{f)df^ it is clear that 
the observed giant increase of Sa/ {a) for ris < (Fig. ^) 
reflects a sudden and dramatic slowing down of the elec- 
tron dynamics. This is attributed to the freezing of the 
electron glass. We also find that, for ris < rig, ^/(/) 
creases exponentially with decreasing T. The observed 
T-dependence of noise (obvious from Fig. ^ is consistent 
with early studies on Si MOSFETs |2^, and it shows 
that the noise in our system cannot be explained by the 
models of thermally activated charge trapping |2^ , ^ , p9| , 
noise generated by fluctuations of T |3^] , noise in the 
hopping regim e M , and in the vicinity of the Ander- 
son transition |l24|. On the other hand, similar increase 
of noise at low T has been observed in mesoscopic spin 
glasses [ ^2|j3^] , in wires in the quantum Hall regime for 
tunneling through localized states Q, and in Si quan- 
tum dots in the Coulomb blockade regime [3^] . 

Another remarkable result, shown in Fig. ||(b), is a 



sharp jump of the exponent a at ris 



9' 



While a ^ 1 



for ris > rig, a ^ l.S below rig, reflecting a sudden shift 
of the spectral weight towards lower frequencies. Simi- 
lar large values of a have been observed in spin glasses 
above the MIT |3^] , and in submicron wires . In gen- 
eral, such noise with spectra closer to 1//^ than to 1// is 
typical of a system far from equilibrium, in which a step 
does not lead to a probable return step. We also have 
the analysis of higher order statistics (non-Gaussianity or 
second spectra [|l2|,|6|) of the noise, showing an abrupt 
change to the sort of statistics characteristic of compli- 
cated multi- state systems just at the density rig at which 
a jumps. This will be described in detail elsewhere. 
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FIG. 5. (a) The averaged noise power spectra Si oc 1/ vs. f for several n^. Solid lines are linear least-squares fits with 
the slopes equal to a. (b) a vs. Us. (c) Relaxation of a following a slow change in Us (lO^^cm"^) from 15.93 to 6.01, carried 
out at T = 0.8 and T = 0.13 K over a period of ~ 4.5 hours each. After a reached a stationary value (A) at 0.8 K, the sample 
was cooled down to 0.13 K. The resulting a (B) differs from a (C) obtained using a different cooling procedure by a factor of 
two, clearly demonstrating history dependent, i. e. non-ergodic behavior. Such behavior is not observed for Us > rig. 



We have demonstrated that the transition to a glassy 
phase is characterized by a sudden, enormous increase 
in the low- frequency noise in a, a sudden shift of the 
spectral weight towards lower /, and a dramatic increase 
of noise with decreasing T. Similar behavior in spin 
glasses was attributed to spin glass freezing [^,^. In 
addition, for < n^, we have observed long relaxation 
times following a large change in Vg^ and history depen- 
dent behavior characteristic of a glassy phase [Fig. ||(c)]. 
We note that in order to obtain reproducible values of 
{a{ns,T)) shown in Fig. ||, it was necessary to vary Us in 
small steps of 4.3 x IQi^cm-^ at the highest T (0.8 K). 

In summary, we present the first evidence of electron 
glass freezing at a well-defined density Ug in a 2D electron 
system in silicon, in agreement with theoretical predic- 
tions [||J^]. Glassy freezing occurs in the regime of very 
low (cr), apparently as a precursor to the MIT. The exis- 
tence of such an intermediate {ric < Us < Ug) glass phase 
is consistent with theoretical predictions pQ| . 
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